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ADD COMMAND TO NETUST TO RUN A SIMULA TION 
EACH TIME A PARAMETER IS ALTERED (i.e. A PLOT) 
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ADD COMMAND TO NETUST TO CALCULATE 
SENSITIVITY OF EACH VECTOR PER EACH PLOT 



ADD COMMAND TO NETUST TO PERFORM 
ADDITIONAL ANA YSIS (OPTIONAL) 
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SENS, sensitivity analysis Simulation Template With Comments: 

^Instruct the netlist builder to show tolerances' 
#tolerance 



It IT 1111 V/VUI 



^Suppress automatic vector saves L ^ 
#nosave \ — 



FIG. 



^Suppress IsSpice4 printout 
#noprint 



\-406 



*Save vectors needed for measurements 
#vector 



\-408 



*Set the output file pointer to the beginning 
* the input net list 
set rewind 



to remove I 

[-410 



*Set the noecho environment for print formatting 
set noecho 



\-412 



*Run the specified simulation and save the results 
#simulation 
set printmode = save 
#mprint 



-414 



^Rename the simulation plot 
nameplot ref 



\-416 



*Set the print format 

SETCOLWK>TH=22 

SETSPICEDIGITS=5 



*Tell the user where we are 

printstatus -t ######## jensitivity Jor_each jarameter_###» 



\-420 



*Loop through all of the parameters 
nextparam null 
while param o null 

* Alter each parameter 

alterparam tolerance(param)/3 

^Simulate 

#simulation 

Save the parameter reference in the new plot 

paranwec = param 

*Tell the user where we are 

printstatus -p paranwec 

*Save the data 
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#mprint 

*Loop through the vectors 
nv = nextvector(null) 
while nvonull 
*Save the sensitivities of scalar data 
if length(nv)= 1 
nv = nv - ref.nv 
end; end if 
*Getthe next vector 
nv = nextvector(nv) 
end; end vector loop 
^Restore the parameter 
unalterparam 
*Get the next parameter 
nextparam 
end; end parameter loop 
*Set print mode for printing output data 
set printmode = print 
unset noecho 
*Loop through the plots 
pi = nextplot(null) 
while pi o null 
if sameplot(ref.default) = 0 
*Loop through the vectors 
nv = nextvector(null) 
while nv o null 
if lengtb(nv) = 1 
*Save the sensitivities 
ifref.nvoO 

nv = (300*NV)/REF.NV 
else 

nv = 3*NV 
end ; end if 
end ; end if 
nv = nextvector(nv) 
end ; end vector loop 
end; end if 
*Get next plot 
pi = nextplot(pl) 
end; end plot loop 
*Tett the user what's going on 
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FIG. 4-2 
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424 



printstatus -t ########sorting sensitivity _for_each _parameter_## 
printstatus -t 

*Make ref the current plot 
setplot ref 

*Sort by descending value 

sort -vd 
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*Loop through the plots 
pi = nextplot(null) 
while pi o null 
if sameplot(ref.default) = 0 

*Print Headers 

SETPARAM PARAMVEC 

printstatus -p paramvec 

ECHO 

ECHO -u M **********SENSmVITY DATA**********" 

ECHO 

ECHO -un "PARAMETER NAME: " 
PRINTNAME PARAMVEC 
ECHO 

ECHO -un " NOMINAL VALUE: " 
PRINTVAL PARAMVEC 
ECHO 
ECHO 

PRINTTEXT -u VECTOR SENSITIVITY% 
ECHO 
ECHO 

*Sort by descending data value 
sort -vd 

*Loop through the vectors and print data 
nv = nextvector(null) 
while nv o null 
if length(nv)= 1 
if ref.nv o 0 
PRINTNAME NV 
PRINTVAL NV 
else 

PRINTNAME NV 
PRINTVAL NV 
ECHO-n* 
end ; end if 
ECHO 
end ; end if 
*Get next vector 
nv = nextvector(nv) 
end; end vector loop 
end ; end if 
*Get next plot 
pl = nextplot(pl) 

end; end plot loop 

ECHO — 
ECHO 

*Print data in output file for SpiceNet to read 
setplot ref 

echo //////////////////// SENSITIVITY analysis Results llliltiltllllllltlllliilll 
#mprint 
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FIG. 4-3 
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C:\spice8d\Circuits\Bandpass.cir Setupl 

*#save V(2) V(3) @C3[i] @C3[p] V(l) @Rl[i] @Rl[p] V(4) 



*#save @C4[i] @C4[p] @R2[i] @R2[p] @R5[i] @R5[p] V(6) V(5) 
*#save @Vl[i] @Vl[p] @V2[i] @V2[p] @V3[i] @V3[p] 

[— 504 



— 502 



*Mias v_4 v(4) 

*#view tran v_4 

TRAN ,05ms 20ms 
.PRINT AC VDB(4) 
.OPTIONS vscale=4 
.PRINT TRAN V 4 



— 506 



FIG. 5B-1 
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C323 10N 
Rl 1 2 10K 
C42410N 
R2201MEG 
R534500K 

X1 03654LM324M{ } 

.SUBCKTLM324M 1 2 345 
* 

CI 11 12 3.000E-12 
C2 6 76.000E-12 
CEE 1099315.8E-15 
DC 5 53DX 
DE 54 5DX 
DLP 9091 DX 
DLN 9290DX 
DP 4 3DX 

EGND99 0POLY(2)3 0 40 0.5.5 

FB 799POLY(5)VBVCVEVLPVLN053.05E6 

+ -50E650E650E6-50E6 

GA 6 0 11 12 37.70E-6 

GCM 0 6 1099 11.92E-9 

IEE 310DC2.476E-6 

HLIM90 0VLIM1K 

Ql 11 213QX 

Q2 12 1 14 QX 

R2 6 9 100.0E3 

RC1 4 11 26.53E3 

RC2 4 1226.53E3 

RE1 13 104.820E3 

RE2 14104.820E3 

REE 109980.78E6 

ROl 8 5 50 

R02 7 99 50 

RP 3 434.71E3 

VB 9 ODCO 

VC 3 53DC2 

VE 54 4DC5.000E-3 

VLIM 7 8DC0 

VLP910DC40 

VLN 0 92 DC 40 

.MODEL DX D(IS=800.0E-18) 

.MODEL QX PNP(IS=800.0E-18 BF=31.58) 

.ENDS 

VI 10 AC=1 PULSE 0-1 IMS 

V250DC=-5 

V3 60DC=5 

.END 
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C:\spice8d\Circuits\Bandpass.cir Setupl 
.OPTIONS vscale=4 
.control 
alias v_4 v(4) 

viewtranv_4 FIG* 6~1 

save v(4) 608 



set rewind 610 

setnoecho 612 

TRAN ,05ms 20ms ~ 

set printmode = save 

echo TRAN Analysis Measurements 

echo 

echo Test 1 Mean 
homeCursors 

print Mean(V(4)) 



614 



nameplotref — 616 

SETCOLWIDTH=22 
SETSPICEDIGITS=5 




618 



printstatus -t 
nextparam null 
while paramomill 

alterparam tolerance(param)/3 

TRAN .05ms 20ms 

paramvec = param 
printstatus -p paramvec 
echo TRAN Analysis Measurements 
echo 

echo Test 1 Mean 
homeCursors 
print Mean(V(4)) 

nv = nextvector(null) 
while nv o null 
iflength(nv) = l 

nv = nv - ref.nv 

end 

nv = nextvector(nv) 
end 

unalterparam 
nextparam 
end 



■jensiti\dtyJor_eachj)arameter_####### — 620 
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set printraode = print 
unset noecho 

pi = nextplot(null) 
while pi o null 
if sameplot(ref.default) = 0 

nv = nextvector(null) 

while nv o null 
if length(nv) = 1 
if refnv o 0 

nv = (300*NV)/REF.NV 
else 

nv = 3*NV 

end 

end 

nv = nextvector(nv) 

end 
end 

pi = nextplot(pl) 
end 



printstatus -t ########sorting jensitivity _for _each _parameter_######## 
printstatus -t 

setplot ref 
sort -vd 

pi = nextplot(null) 
while pi o null 
if sameplot(ref.default) = 0 

SETPARAM PARAMVEC 

printstatus -p paramvec 

ECHO 

ECHO -u "**********SENSITIVITY DATA**********" 

ECHO 

ECHO -un "PARAMETER NAME: " 
PRINTNAME PARAMVEC 
ECHO 

ECHO -un H NOMINAL VALUE:" 
PRINTVAL PARAMVEC 
ECHO 
ECHO 



FIG. 
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PRINTTEXT -u VECTOR SENSITMTY% 

ECHO 

ECHO 

sort -vd 

nv = nextvector(nulI) 

while nv o null 
if length(nv) = 1 
ifref.nvoO 

PRINTNAMENV 
PRINTVALNV 
else 

PRINTNAMENV 

PRINTVALNV 

ECHO-n* 

end 
ECHO 

end 

nv = nextvector(nv) 

end 



end 

pi = nextplot(pl) 
end 

ECHO 
ECHO 
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FIG. 



628 



setplotref 

echo //////////////////// SENSITIVITY analysis Results ////////////////////////// 

echo TRAN Analysis Measurements 

echo 

echo Test 1 Mean 
homeCursors 

print Mean(V(4)) 



— 630 



.endc 

C323 10NTOL=5% 

R11210KTOL=2% 

C42410NTOL=5% 

R2201MEGTOL=2% 

R534500KTOL=2% 

X1 03 654LM324M{ } 
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.SUBCKTLM324M 1 2345 
* 

CI 11 12 3.000E-12 
C2 6 7 6.000E-12 
CEE 1099315.8E-15 
DC 553DX 
DE 54 5DX 
DLP 9091 DX 
DLN 9290DX 
DP 4 3DX 

EGND99 0POLY(2)3 0 40 0.5.5 

FB 799POLY(5)VBVCVEVLPVLN053.05E6 

+ -50E650E650E6-50E6 

GA 6 011 1237.70E-6 

GCM 0 6 1099 11.92E-9 

IEE 310DC2.476E-6 

HLIM90 0VLIM1K 

Ql 11 213QX 

Q2 12 114QX 

R2 6 9100.0E3 

RC1 4 1 1 26.53E3 

RC2 4 1226.53E3 

RE1 13 104.820E3 

RE2 14 104.820E3 

REE 109980.78E6 

ROl 8 5 50 

R02 79950 

RP 3 434.71E3 

VB 9 0DC0 

VC 353DC2 

VE 54 4DC5.000E-3 

VLIM 7 8DC0 

VLP 91 ODC40 

VLN 092 DC 40 

.MODEL DXD(IS=800.0E-18) 

.MODEL QX PNP(IS=800.0E-18 BF=31.58) 

.ENDS 

VI 1 0AC=1 PULSE 0-1 IMS 

V250DC=-5 

V360DC=5 

.END 
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**********SENSITIV1TY DATA*****"*** 

PARAMETER NAME: r5 
NOMINAL VALUE: 500.00K 

VECTOR SENSITIV]TY% 

mean(v(4)) 1.5111 

**********SENSITIVTTY DATA**"****** 

PARAMETER NAME: r2 
NOMINAL VALUE: l.OOOOMeg 

VECTOR SENSITIVITY 0 /.. 

mean(v(4)) 17.265M 

**********SENSmVITYDATA********** 

PARAMETER NAME: c4 
NOMINAL VALUE: 10.0000N 

VECTOR SENSITIV1TY% 

mean(v(4)) -752.77M 

**********SENSnTVlTY DATA**"****** 

PARAMETER NAME: rl 
NOMINAL VALUE: 10.0000K 

VECTOR SENSITIVITY% 

mean(v(4)) -571.46M 

**********SENSITIVITY DATA********** 

PARAMETER NAME: c3 
NOMINAL VALUE: 10.0000N 

VECTOR SENSITIVITY% 



mean(v(4)) 4.5201 



sensitivity analysis results 
tran analysis measurements 

test 1 mean 

mean(v(4)) = 2.086052e-001 
Total run time: 0.583 seconds. 
Total run time: 0.583 seconds. 



FIG. 7 



Memory remaining = 1996210 Kbytes 
Memory Used = 14401 Kbytes 
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RSS, root summed square analysis Simulation Template With Comments: 



instruct the netlist builder to show tolerances 
#tolerance 



*Suppress automatic vector saves 
#nosave 



FIG. 8-1 



^Suppress IsSpice4 printout \-ggg 
#noprint | 



*Save vectors needed for measurements 
#vector 



*Set the output file pointer to the beginning to remove J 

* the input net list 1-810 

set rewind | 

*Setthe noecho environment for print formatting I gj2 
set noecho I — 



*Run the specified simulation and save the results 

#simulation 

set printmode = save 

#mprint 



} 



814 



*Set the print format 
SETCOLWIDTH=22 V-818 
SETSPICEDIGITS-5 _| " 



^Rename the simulation plot ~~ I ^ 
nameplot ref | — 



*Loop through all of the parameters 

nextparam null 

*Tell the user where we are 

printstatus -t " //////////////// sensitivity for each parameter 

while paramo null 

% Alter each parameter 

alterparam tolerance(param)/3 

^Simulate 

#simulation 

*Save the parameter reference in the new plot 

paramvec = param 

*Tell the user where we are 

printstatus -p paramvec 

*Save the data 

#mprint 



820 



-822 
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*Loop through the vectors 
nv = nextvector(null) 
while nv o null 

*Save the sensitivities of scalar data 

iflength(nv) = 1 
nv = nv-ref.nv 

end; end if 

*Get the next vector 

nv = nextvector(nv) J- 822 

end; end vector loop 

^Restore the parameter 

unalterparam 

*Getthe next parameter 

nextparam 
end ; end parameter loop 
*&/ print mode for printing output data 
set printmode = print 
unset noecho 

*Make a new plot for rss data and make it the current plot 

newplot rss ref.default ref.default 

setplotrss 

*Loop through the vectors 
nv = nextvector(null) 
while nv o null 

^Initialize the scalar data to zero 

IF LENGTH(NV) = 1 
nv = 0 

END;endif 

* Get the next vector 

nv = nextvectoifnv) 

end; end vector loop 

*Loop through the plots 

pl = nextplot(null) 

*Tell the user where we are 

printstatus -t " //////////////// computing rss for each measurement 



FIG. 8-2 



~^824 



-832 



-826' 



while pi o null 
*Select plots 

if sameplot(rss.default) = 0 
if sameplot(ref.default) = 0 

*Tett the user what we are doing 
printstatus -1 pi 

*Mahe the saved parameter reference, paramvec current 
SETPARAM PARAMVEC 
*Get the next vector 
nv = nextvector(null) 
*Print formatted data 
ECHO 
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ECHO -U "****»***** RSS DATA *********»" 

ECHO 

ECHO -UN " PARAMETER NAME: " 
PRINTNAME PARAMVEC 
ECHO 

ECHO -UN " NOMINAL VALUE: " 
PRINTVAL PARAMVEC 
ECHO 

ECHO -UN "TOLERANCE VALUE: " 
PRINTTOL PARAMVEC 
ECHO 
ECHO 

PRINTTEXT -UN VECTOR 

PRINTTEXT -U SENSITIV1TY% RSS_CONTRIBUTION 

ECHO 

ECHO 

*Loop through vectors 
while nv o null 
if length(nv) = 1 
IFREF.NVO0 
*Calculate the RSS percentage if value is not zero 
PRINTNAME NV 
NEWNV = (300*NV)/REF.NV 
PRINTVAL NEWNV 
ELSE 

*Calculate the RSS if value is not zero 
PRINTNAME NV 
NEWNV = 3*NV 
PRINTVAL NEWNV 
ECHO-N* 
END;endif 

*Save and print the sum of squares 
rss.nv = rss.nv + nv * nv 
PRINTVAL RSS.NV 
ECHO 
end ;end if 
* Get next vector 
nv=nextvector(nv) 
end ;end vector loop 
* Sort plot by descending value 
sort-vd 
end ;end if 
end ;end if 
* Get next plot 
pl = nextplot(pl) 

end ;end plot loop 
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FIG. 8-3 
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*Sort the rss plot by descending value 
setplotrss 
SORT-VD 
*Print Headers 
ECHO 

ECHO -U "«*********RSS H1/L0 ANALYSIS RESULTS**********" 

ECHO 

PRINTTEXT -UN VECTOR 
SETCOLWK)TH=15 

PRINTTEXT -U NOMINAL RSS-VALUE TOLERANCE% HI_VALUE L0_VALUE 

ECHO 

ECHO 

*Make a new plot for results 
newplot hires ref.default ref.default 
*Loop through the vectors 
nv = nextvector(null) 
while nv o null 
if length(nv) = 1 

*Print formatted data 
SETCOLWIDTH=22 

PRINTNAMENV 
SETC0LWIDTH=15 
PPJNTVALREF.NV 

NV = 3*SQRT(ABS(NV)) 
PRINTVALNV 
IF REF.NV o 0 
NEWNV = (100*NV)/REF.NV 
ELSE 

NEWNV = NV*0 
END 

PRINTVAL NEWNV 
NV = REF.NV + NV 

HI_RSS = REF.NV + NV 

L0_RSS = REF.NV-NV 

PRINTVAL HI_RSS 

PRINTVAL L0_RSS 

ECHO 
end ; end if 
*Get next vector 
nv = nextvector(nv) 

end; end vector loop 

ECHO 
ECHO 

*Print data in output file for SpiceNet to read 
setplotrss 

echo llllllilllllllllilll RSS HI analysis Results i 
#mprint 

RUSAGE ELAPSED 



834 



\-m FIG. 8-4 
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EVA, Extreme Value Analysis Simulation Template With Comments: 

^Instruct the netlist builder to show tolerances ' 



#tolerance 

^Suppress automatic vector saves 
#nosave 



[-902 



~^-904 



* Suppress IsSpice4 printout I ^ 
#noprint | — 



*Save vectors needed for measurements 
#vector 



FIG. 9-1 



Set the noecho environment for print formatting 

set rewind 910 

set noecho 912 



*Run the specified simulation and save the results ' 
#simulation 

pltype = 0 jdentify the plot type for later use 

set printmode = save 

#mprint 



-914' 
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*Set the print format 
SETCOLWIDTH=22 
SETSPICEDIGITS=5 
nameplotref — 916 
newplot evahi ref.default refdefault 
evahi.pltype = 0 Identify the plot type for later use 
*Print status for the user 

printstatus -t "######## sensitivity for each parameter 
*Loop through the parameters 
nextparam null 
while param o null 

% Alter each parameter 

alterparam tolerance(param)/3 

^Simulate, making a new plot for results 

#simulation 

*Save the current parameter reference 

paramvec = param 

*Tell the user where we are 

printstatus -p paramvec 

*Save the tol and paramval 

paramtol = tolerance(param) 

paramval = getparam(param) 

pltype = 1 Identify the plot type for later use 

*Save the simulation results 



~^-932 



7^-920 
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#mprint 

*Loop through all the vectors 
nv = nextvector(null) 
while nv o null 

*Save the sensitivities for all scalar measurements except pltype 

if length(nv) = 1 
if nv o pltype 

nv = nv - ref.nv 
end ;end if 

end ;end if 

nv = nextvector(nv) 
end ;end vector loop 
^Restore the param 
unalterparam 
*Getthe next param 
nextparam 

endjend parameter loop 

*Makeref the current plot 
setplot ref 

*Loop through the vectors in ref 
nv = nextvector(null) 
*Tell the user where we are 
printstatus -t " Mi II Ml Ml II measurements M M Ml M l I I " 
while nv o null 
iflength(nv) = 1 

*Loop through all the plots containing scalar vectors 

pi = nextplot(null) 

while pi o null 
if pltype = 1 

* the inner loop, we are looping through each sensitivity plot looking at the same vector 

* we will alter the parameter id'd by paramvec to maximize/minimize the vector 

setparam paramvec 

^Change each parameter to its worst case extreme value 
ifnv >=0 

alterparam paramtol 
else 

alterparam -paramtol 
end 
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pi = nextplot(pl) 
end 

^Simulate for the extreme case and save the data in a new plot 

#simulation 

#mprint 

pltype = 2 identify the plot type for later use 

* if we want sensitivity at the extreme, we need to go through each param 

* and change it to be a bit different than it is at the extreme, run a simulation , 
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* and mark it as pltype 3 along with its paramvec, then we can take the diff 
*from the pltype = 2 to get the sensitivity at the extreme if the sign at the extreme 

* is different than at the nominal we can report an error or go on to do worst case 
*for worst case, we need to reduce the param change by 1/2 and do this over again. 

* either continue in this loop or make a wc loop afterward., save the paramvalue 

* and tolerance 

evahlnv = nv 
ifnvo pltype 
*Tell the user where we are 
printstatus -n nv 
end ;endif 
end ;end plot loop 
; get the next vector 
nv = nextvector(nv) 

end ;end vector loop 

*Set print mode and print header 
set printmode = print 
unset noecho 
setplotEVAHI 
ECHO 

ECHO -U ''**********EVA PARAMETER LIST**********" 
ECHO 

PRINTTEXT -UN PARAMETER 
PRINTTEXT -U NOMINAL TOLERANCE 
ECHO 
ECHO 

*Loop through the parameters 
nextparam null 
while param o null 

*Extract the saved param reference and print its data 

paramvec = param 

PRINTNAME PARAMVEC 

PRINTVAL PARAMVEC 

PRINTTOL PARAMVEC 

ECHO 

*getthe nextparam 
nextparam 
end ;end param loop 

ECHO iprint a blank line 

*Make a new plot to hold sorted results 
newplot evasort ref.default refdefault 
*Make ref the current plot 
setplot REF 

*Loop through all vectors in ref 
nv = nextvectorfnull) 
while nv o null 
if length(nv) = 1 j 
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*save the result in evasort as a percent of its value J 
if refnv o 0 

evasortnv = ((evahi.nv-ref.nv)*100)/ref.nv 
else 

evasort.nv = 0; 
end ;end if 
end ;end if 
*Get the next vector 
nv = nextvector(nv) 

end ;en d vector loop 

*Print some headers 
ECHO 

ECHO -U "**********EVA-HI RESULTS**********" 

ECHO 

PRINTTEXT -UN VECTOR 

PRINTTEXT -U NOMINAL EVA-HI CHANGE% 

ECHO 

ECHO 

setplot evasort 

*Sort evasort by descanting data 
sort-VD 

*Loop through the vectors 
nv = nextvector(null) 
while nv o null 
*Ifits the correct data in the correct plot, print it 
if length(nv) = 1 
ifnvopltype 
PRINTNAMENV 
PRINTVAL REF.NV 
PRINTVAL EVAHLNV 
PRINTVAL EVASORT.NV 
ECHO 
end ; end if 
end ; end if 
*Get the next vector 
nv = nextvector(nv) 
end ; end vector loop 
ECHO 

ECHO 



* now the eva results are in pltype = 2 plots 

*Print the results so SpiceNet can read the eva-hi data 

set printmode = print 

unset noecho 

setplot evahi 

echo liiilllllllllllllll! EVA HI analysis Results ////////////////////////// 
#mprint 
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WCS, Worst Case by Sensitivity Simulation Template With Comments: 

^Instruct the netlist builder to show tolerances j_ 
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#tolerance 



^Suppress automatic vector saves I 
#nosave I 



FIG. 10-1 



^Suppress IsSpice4 printout L 
#noprint — I 
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*Save vectors needed for measurements I 
#vector | 

*Set the output file pointer to the beginning to remove 
* the input net list 
set rewind 



*Setthe noecho environment for print formatting I 

set noecho I 



*Run the specified simulation and save the results 

#simulation 

set printmode = save 

#mprint 
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*Set the print format 

SETCOLWIDTH=22 

SETSPICEDIGITS=5 
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^Rename the simulation plot I jq j$ 
nameplotref I 



*Make a newplot for results 
newplot result ref.default ref default 
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*Set the current plot to ref — 1 ^4 
setplotref _| 

*P rint status for the user 

printstatus -t "######## sensitivity for each parameter 



*Loop through the parameters 

nextparam null 

while paramo null 
* Alter each parameter 
alterparam tolerance(param)/3 
^Simulate, making a new plot for results 

#simulation 
*Save the current parameter reference 
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pararavec = param 
Vnform the user about what's being done 

printstatus -p paramvec 
*Make and save the measurements 
#mprint 

*Loop through the vectors 
nv = nextvector(null) 
while nv o null 

*Save the sensitivity of scalar quantities 
iflength(nv)= 1 
nv = nv - refnv 
*Save the worst case -hi value 
resultnv = resultnv + abs(3*nv) 
end ;end if 
nv = nextvector(nv) 
end ;end vector loop 
*restore the parameter value 
unalterparam 
*get the next parameter 
nextparam 

end ;end parameter loop 
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FIG. 10-2 
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*Set the print mode to print instead of save 

set printmode = print 

^Restore the echo mode for printing 

unset noecho 

*Set result to the current plot 
setplot result 
*Print the header 
ECHO 

ECHO -U ''**********WCS PARAMETER LIST**********' 
ECHO 

PRINTTEXT -UN PARAMETER 
PRINTTEXT -U NOMINAL TOLERANCE 
ECHO 
ECHO 

*Loop through the parameters 
nextparam null 
while param o null 
paramvec = param 
*Print the row 

PRINTNAME PARAMVEC 

PRINTVAL PARAMVEC 

PRINTTOL PARAMVEC 
ECHO 
nextparam 
end 

ECHO 
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*Make a new plot to hold sorted results I 25/25 

newplot wcsort ref.default ref.default 
*Set the current plot to ref 
setplot REF 

*Loop through its vectors 
nv = nextvector(null) 
while nv o null 
^Calculate the wcasa percent change results 
iflength(nv)= 1 
ifref.nvoO 

wcsortnv = ((result.nv-ref.nv)*100)/ref.nv 
else 

wcsortnv = 0; 
end ;end if 
end ;end if 
nv = nextvector(nv) 

end ;end vector loop 

*Print headers 

ECHO 

ECHO -U "**********WCS-HI RESULTS**********" 

ECHO 

PRINTTEXT -UN VECTOR 
PRINTTEXT -U NOMINAL WCS-HI CHANGE% 
ECHO 
ECHO 

;sort wcsort by descending value 
setplot wcsort 
sort-VD 

*Print the ordered list 
nv = nextvector(null) 
while nv o null 
iflength(nv)= 1 

PRINTNAMENV 

PRINTVAL REF.NV 

PRINTVAL RESULT.NV 

PRINTVAL WCSORT.NV 

ECHO 
end 

nv = nextvector(nv) 

end 
ECHO 

ECHO 

*Set the current plot to the wc results 
setplot result 

echo //////////////////// WCS HI analysis Results ////////////////////////// 
*Print the measured results in a form that can be read back into SpiceNet |— 1030' 
#mprint 

*R eport the elapsed time in the output file 

rusage elapsed 



- 1040 



